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1.  INTRODUCTION 


A  bibliography  (Sherman  1981)  has  recently  been  prepared  listing  work  related  to  the 
occurrence  of  gusts  of  aeronautical  significance  in  Australia.  Only  a  small  number  of  investi¬ 
gations  give  a  decomposition  by  altitude  of  turbulence  encountered  in  routine  flying  operations. 
These  investigations  are  listed  in  Table  1,  and  the  distances  flown  in  each  height  band  are 
indicated  graphically,  in  Figure  1,  by  rectangles  whose  area  is  proportional  to  the  distance 
flown  in  the  height  band  multiplied  by  a  weighting  factor  which  allows  for  the  relevance  of  the 
flying  to  general  Australian  aviation.  The  weighting  factors  are  given  in  column  6  of  Table  1 . 

The  data  for  three  of  the  programs  (NZ  Viscount,  Comet  and  707)  are  either  not  available 
in  the  published  literature  or  have  not  been  published  separately  for  the  Australian  region.  The 
data  are,  however,  stored  on  magnetic  tape  at  the  RAE  and  could  probably  be  obtained  in  a 
suitable  form.  However,  Kaynes  has  indicated,  in  response  to  a  request  for  the  data,  that  the 
age  of  the  magnetic  tapes  and  changes  in  the  RAE  Computer  System  present  potential  problems 
for  extracting  valid  data.  The  data  from  the  Mirage  V-g-h  program  are  still  being  analysed,  but 
they  are  probably  of  dubious  general  application  because  of  the  special  highly  responsive 
characteristics  of  fighter  aircraft. 

Considerably  more  V-g-h  data  was  acquired  from  the  Boeing  727  aircraft  than  the  600 
hours  analysed  by  Hunter  and  Fetner  (1967).  In  fact,  at  DCA’s  request,  NASA  fitted  V-g-h 
recorders  into  a  727  from  each  fleet  (Ansett  and  TAA)  from  about  1964  to  April  1968,  and 
into  a  DC-9  from  each  fleet  from  about  1966  to  1970.  In  a  paper  discussing  “initial”  data. 
Hunter  (1967)  indicated  that  533  hours  of  Ansett  727  data  and  346  hours  of  TAA  727  data  had 
been  analysed  at  that  stage.  However,  the  only  altitude  decomposition  available  is  in  Hunter 
and  Fetner’s  (1967)  paper  on  the  first  600  hours.  No  results  are  available  from  the  DC-9  program. 
The  complete  results  of  these  two  programs  would  be  a  very  significant  addition  to  the  total 
knowledge  of  the  Australian  environment,  but  the  analysis  does  not  appear  to  have  been  pur¬ 
sued  at  the  time,  and  as  NASA  have  now  closed  their  V-g-h  program  it  is  not  certain  whether 
the  Australian  Department  of  Transport  (which  now  includes  the  former  DCA)  would  be  able 
to  obtain  even  the  raw  data. 

The  results  of  the  programs  for  which  data  are  available  have  been  summarised  in  the 
Appendix.  In  all  this  data  there  are  really  only  three  major  programs  covering  the  Australian 
area.  These  are  the  Qantas  Super  Constellation  flights,  mainly  around  15,000  ft,  the  Viscount 
flights,  mainly  around  20,000  ft,  and  the  Boeing  727  flights,  which  were  mainly  around 
35,000  ft. 


2.  THE  VARIATION  OF  TURBULENCE  WITH  ALTITUDE 

One  of  the  well  recognised  design  codes  for  gust  loads  has  been  produced  by  the  Royal 
Aeronautical  Society  of  London  as  the  ESDU  data  item  69023.  We  will  compare  the  available 
Australian  design  data  with  the  design  curves  shown  in  that  data  item.  The  distance  between 
gusts  which  exceed  3  m/s  is  conventionally  denoted  tio  * 

Figures  2  and  3  show  the  variation  of  t io  with  altitude,  for  the  data  listed  in  Table  I,  and 
on  the  same  graphs  the  design  curves  from  Figures  1  and  2  of  the  ESDU  data  item  are  shown. 
(In  these  and  the  later  figures  the  octagons  around  each  data  point  have  8n  area  proportional 


*  Following  present  aeronautical  usage  we  have  retained  the  imperial  units  of  nautical 
miles  for  horizontal  distance  and  feet  for  altitude.  Gust  velocities  have  been  converted  from 
ft/s  to  m/s  using  the  approximate  conversion  factor  0-3  ms  =  1  ft/s.  Following  the  ESDU 
data  item  we  have  retained  the  nomenclature  ((to  rather  than  (s)  for  the  distance  between  gusts 
which  exceed  3  m/s  (i.e.  10  ft/s). 
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to  the  distance  flown  in  the  altitude  band,  multiplied  by  the  appropriate  weighting  factor  from 
column  6  of  Table  I .)  At  low  attitudes  (below  5000  ft)  there  is  insufficient  data  to  make  any 
conclusive  statement.  The  Australian  Bristol  Freighter  data  indicate  a  low  incidence  of 
turbulence,  which  might  be  expected  since  much  of  the  flying  was  over  the  sea  (across  Bass 
Strait).  Likewise  the  New  Zealand  Bristol  Freighter  data  for  the  low  altitudes  indicate  low 
turbulence,  and  much  of  this  data  represents  flying  across  Cook  Strait.  The  medium  altitude 
turbuience  tends  to  average  out  near  the  design  curve;  the  Viscount  data  being  milder  than 
wouid  be  expected,  given  the  fact  that  the  aircraft  were  not  equipped  with  radar,  and  the  Super 
Constellation  data  for  the  Australian  region  being  more  severe  than  would  be  expected.  At 
high  altitudes  (around  35,000  ft)  the  TAA  727  data  is  a  little  less  severe  than  the  appropriate 
design  curve  (for  aircraft  with  cloud  warning  radar)  whilst  the  Ansett  data  is  much  more  severe 
than  even  the  design  curve  for  aircraft  without  radar.  A  relatively  severe  turbulence  environ¬ 
ment  at  high  altitude  is  quite  credible  because  this  element  of  Australia’s  flying  occurs  near  the 
lattitude  and  altitude  of  the  jet  stream  core,  and  predominantly  over  mountainous  terrain. 
Therefore,  it  seems  to  me  highly  desirable  to  obtain  more  data  for  this  altitude  band. 


3.  THE  RATIO  OF  UPGUSTS  TO  DOWNGUSTS 

Following  Bullen  (1966)  the  ESDU  data  item  69023  has  indicated  that  at  low  altitudes 
positive  acceleration  increments  occur  more  frequently  than  negative  increments.  This  is 
frequently  attributed  to  the  effect  of  manoeuvre  loads  (see,  for  example,  Bullen,  op.  cit.).  The 
ratio  of  upgusts  to  downgusts  for  the  Australian  data  has  been  plotted  in  Figure  4,  which 
parallels  Figure  4  of  the  ESDU  data  item.  The  curve  shown  is  computed  by  the  expression 

0  85//  +  9100 
//  +  3300 

for  altitudes,  H  (ft),  less  than  38,667  ft,  and  is  taken  as  unity  for  greater  altitudes.  This  curve 
is  a  close,  but  not  exact,  tit  to  the  one  shown  in  the  ESDU  data  item.  It  can  hardiy  be  said  that 
the  data  support  the  curve,  but  neither  do  they  suggest  any  other  curve. 


4.  THE  VARIATION  OF  GUST  FREQUENCY  WITH  GUST  MAGNITUDE 

Figures  5  and  6  show,  for  three  different  altitude  bands,  the  frequency  of  occurrence  of  gusts 
of  various  magnitudes,  relative  to  the  number  of  times  a  3  m/s  gust  is  exceeded.  The  Figures 
parallel  Figures  5  and  6  in  the  ESDU  data  item,  and  include  the  same  design  curves  as  the  ESDU 
data  item.  Figure  5  shews  the  design  curves  for  an  aircraft  not  equipped  with  cloud  warning 
radar  whilst  Figure  6  shows  the  design  curves  for  an  aircraft  with  cloud  warning  radar. 

The  data  for  the  low  altitude  band  (0-5000  ft)  fall  very  close  to  the  ESDU  curve  for  air¬ 
craft  without  cloud  warning  radar.  The  greater  part  of  these  data  was  obtained  from  the  Bristol 
v  md  the  TAA  Viscount,  none  of  which  were  equipped  with  radar,  so  it  appears  that 
this  one  of  the  ESDU  design  curves  is  well  substantiated  for  Australian  conditions. 

The  data  for  the  medium  altitude  (15,000-20,000  ft)  band  come  mainly  from  the  Super 
Constellation  and  Viscount,  neither  of  which  were  equipped  with  cloud  warning  radar.  The 
high  magnitude  gusts  appear  to  occur  a  little  less  frequently  relative  to  the  3  m/s  gusts  than 
would  be  expected  from  the  ESDU  data  curves.  Several  explanations  are  possible: 

(o)  Severe  gusts  tend  to  occur  mainly  in  thunderstorms,  which  the  Qantas  Super  Constel¬ 
lation  had  greater  freedom  to  avoid  near  Australia  than  in  the  more  congested  skies  of 
the  northern  hemisphere. 

(b)  Much  of  the  data,  being  from  the  Super  Constellation,  may  have  been  obtained  from 
flights  on  a  regular  schedule  which  -outinely  passed  through  Australian  air  space  at  a 
time  of  day  when  thunderstorm  activity  was  minimal. 

(c)  Much  of  the  Super  Constellation’s  data  were  from  flights  over  the  sea.  The  ESDU 
design  curves  indicate  that  terrain  effects  on  the  incidence  of  3  m/s  gusts  are  negligible 
above  about  19,000  ft.  It  is  possible,  however,  that  the  stronger  gusts  do  show  the 
effect  of  terrain  to  higher  altitudes. 


3 


(<0  Much  of  the  Qantas  flying  it  in  low  latitudes.  There  may  be  a  latitudinal  effect  with 
strong  gusts  occurring  relatively  less  frequently  in  tropical  regions  than  in  mid-latitudes, 
(e)  The  differences  may  only  be  a  result  of  the  small  data  sample.  About  four  times  the 
distance  was  flown  at  mid-altitudes  as  at  low  altitudes,  but  since  the  expected  number 
of  gusts  per  unit  distance  is  only  2%  of  what  is  expected  at  low  altitudes,  the  expected 
number  of  gusts  encountered  at  mid-altitudes  is  only  8%  of  that  for  low  at  t.udes. 

The  data  for  the  high  altitude  (30,000-40,000  ft)  band  are  only  available  up  to  the  6  m/s 
level  because  of  the  relatively  low  incidence  of  turbulence  at  high  altitudes.  Such  data  as  exist 
fail  closer  to  the  severe  (no  radar)  design  curve  than  to  the  milder  (aircraft  equipped  with  radar) 
design  curve.  This  gives  further  weight  to  the  opinion  expressed,  at  the  end  of  Section  2,  that  it 
is  desirable  to  obtain  more  data  for  aircraft  operating  at  high  altitudes. 


5.  CONCLUSIONS 

(1)  Most  of  the  flight  data  available  are  for  scheduled  passenger  transport  operations  for 
which  the  Australia-wide  incidence  of  turbulence  can  at  best  be  defined  to  within  a  factor  of 
about  3,  and  at  many  heights  the  uncertainty  is  much  greater. 

(2)  From  the  Australian  data  it  is  not  possible  to  make  any  distinction  between  types  of 
flying  (e.g.  cruise,  climb  and  descent  or  special  missions)  nor  is  there  much  information  on  the 
higher  or  lower  incidence  of  turbulence  to  be  expected  in  various  regions  of  Australia. 

(3)  The  available  data  are  not  incompatible  with  the  ESDU  data  item  69023,  except  that  at 
high  altitudes  (above  30,000  ft)  there  is  some  indication  that  turbulence  is  encountered  more 
frequently  than  would  be  predicted  by  the  ESDU  design  curves,  and  particularly  so  for  the 
stronger  gusts. 

(4)  A  greater  than  average  frequency  of  high  altitude  turbulence  might  be  expected  because 
the  bulk  of  Australia's  scheduled  flying  is  in  the  Sydncy-Canberra-Melbourne  region  where  there 
is  both  a  strong  jet  stream  and  a  significant  mountain  range.  This  expectation  reinforces  the 
significance  of  the  indications  from  the  available  data  and  so,  in  view  of  the  paucity  of  data 
above  30,000  ft,  it  seems  desirable  to  acquire  more  data  at  high  altitudes. 
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A.O  General 

In  order  to  bring  all  the  data  to  a  common  basis,  various  calculations  have  been  performed 
for  each  separate  program  as  indicated  in  the  later  sections  of  the  Appendix.  A  frequent  problem 
was  to  interpolate  tables  of  exceedances  for  gust  velocities  other  than  the  measured  values.  This 
was  done  by  assuming  that  between  two  measured  values  of  gust  velocity,  the  logarithm  of  the 
number  of  exceedances  varied  linearly  with  the  gust  velocity.  A  small  degree  of  extrapolation 
was  occasionally  performed,  usually  by  assuming  that  in  a  table  of  numbers  of  observed  exceed¬ 
ances,  the  first  zero  value  was  replaced  by  0-1. 

There  arc  certain  observations  which  have  been  made  many  times  before  in  connection 
with  various  programs  of  measurements  of  turbulence  occurrence.  They  are  reproduced  here 
(in  general  without  acknowledgement)  because  they  carry  a  warning  about  how  the  data  may  be 
interpreted. 

(a)  The  occurrence  of  turbulence  is  not  simply  an  overall  (space-time)  atmospheric  average 
because,  in  aircraft  design  or  flight  monitoring,  it  is  the  loads  actually  experienced  by 
the  aircraft  which  matter,  and  these  depend  on  details  of  aircraft  utilization  (which 
may  depend  for  example  on  the  existence  of  weather  radar  or  on  the  turbulence  avoidance 
techniques  used  by  the  pilot).  Several  of  the  following  observations  are  particularizations 
of  this  general  one. 

( b )  Cruise  conditions  are  generally  less  severe  than  climb  and  descent  because,  with  some 
limitations,  a  pilot  can  choose  a  cruise-height  of  minimum  turbulence. 

(c)  Initial  climb  and  final  descent  are  usually  more  severe  than  other  climb  and  descent 
because,  being  generally  in  controlled  air  space,  a  pilot  has  less  freedom  in  turbulence 
avoidance  than  in  other  locations. 

(</)  In  some  cases  it  has  been  suggested  (Bullen  1966)  that  the  decrease  of  gusts  with  altitude 
is  related  to  the  fact  that  pilots  sometimes  climbed  to  higher  altitudes  than  they  would 
otherwise  have  chosen  in  order  to  avoid  turbulence. 

(,  )  Counting  accelerometer  data  have  generally  shown  that  at  low  altitudes  upgusts  occur 
more  frequently  than  down  gusts.  This  may  be  due  to  the  presence  of  manoeuvres  which 
are  not  distinguished  from  turbulence  in  counting  accelerometer  records.  However, 
most  of  the  manoeuvres  made  by  transport  aircraft  appear  sufficiently  mild  to  not 
influence  the  number  of  significant  loads.  It  is  possible  that  the  excess  of  updrafts  is  a 
genuine  atmospheric  phenomenon  because,  in  convective  fields,  small  concentrated 
updrafts  tend  to  occur  in  a  field  of  air  which  elsewhere  is  settling  slowly  and  uniformly. 
There  is,  however,  a  problem  still  requiring  explanation.  Almost  all  of  the  early  data 
(DC3,  DC4  or  Bristol  Freighter)  measured  with  V-g  recorders  in  Australia  (see  Hooke 
1956)  show  downgusts  to  be  very  much  more  frequent  than  upgusts. 

(/)  There  is  one  source  of  difference  between  data  from  US  gust  measurements  with 
with  V-g-h  recorders,  and  UK  measurements  with  counting  accelerometers  which  is 
additional  to  the  other  differences  remarked  on  by  Bullen  (1958).  Counting  accelerometers 
are  unable  to  distinguish  between  gusts  and  manoeuvres,  but  V-g-h  recorders  can,  and 
in  some  cases  (e.g.  Hunter  and  Fetner  1967)  the  effect  of  manoeuvres  has  been  removed 
from  US  gust  data.  With  transport  aircraft  such  as  those  summarised  here,  the  effect 
is  slight,  but  with  smaller  aircraft  the  differences  may  be  considerable. 


{g)  Even  fairly  long  samples  of  data  may  be  influenced  by  regular  scheduling,  especially  if 
certain  routes  are  only  flown  once  per  day.  For  example,  Bacon  (I9S6)  found  a  low 
incidence  of  turbulence  on  the  Adelaide-Perth  route  which  he  attributed  to  the  fact 
that  the  regular  scheduled  flight  left  Adelaide  late  in  the  afternoon  and  arrived  at  Perth 
in  the  evening.  The  return  flight  was  early  the  next  morning  so  several  of  the  low  altitude 
segments  were  flown  near  night  time  hours  when  turbulence  might  be  expected  to  be  low. 

(A)  Some  minor  variations  between  programs  occur  because  in  some  cases  conversion 
between  acceleration  and  derived  gust  velocity  is  based  on  fixed  average  values  of  air¬ 
craft  weight  and  lift-curve  slope,  whilst  other  more  elaborate  analyses  allow  for  variations 
of  these  factors  with  height,  velocity  and  time  during  the  flight. 


A.l  Bristol  Freighter  (ANA) 

Reference:  Baum  and  Hooke  (1933) 

The  reference  gives  a  table  of  occurrences  of  various  acceleration  ranges  in  10  altitude 
bands.  These  were  converted  to  exceedances  of  Un  using  the  following  formula  and  parameter 
values : 

2(W/S) 


Ude 


PoVeOFo 


An 


Dde  =  Derived  equivalent  gust  velocity 
W  —  Aircraft  weight 
S  —  Wing  area  (1487  sq.  ft) 

WIS  =26-8  Ib/sq.  ft 

po  =  Sea-level  air  density  (0-002378  slug/cu.  ft) 

Ve  —  Equivalent  air  speod  (150  kn=253  ft/sec) 
a  =  Lift  curve  slope  (4-4/rad) 

Fo  —  Gust  alleviation  factor  obtained  from  Figure  7  of  ESDU  data  item  69023  using  the 
following  additional  data 
=  Aspect  ratio  (7-8) 

/<  —  Mass  parameter  [2(  W/S)/gpae] 
c  =  Mean  aerodynamic  chord  (13-8  ft) 


For  three  different  heights  we  obtain: 


Height 

(ft) 

p 

(slug/ft3) 

t 

Fo 

Viol  An 
(ft/sec/g) 

0 

0-002378 

US 

30-2 

3,000 

0-002049 

13-4 

28-9 

10,000 

0-001756 

15-6 

mm 

27-7 

For  other  heights  the  value  of  l/<te/An  was  interpolated  linearly.  The  resultant  table  of  gust 
velocity  exceedances  is  shown  as  Table  A1  at  the  end  of  this  Appendix. 


A.2  Viscount  (TAA) 

Reference:  Bacon  (1956) 

The  reference  gives  a  graph  (Fig.  I)  of  exceedances  of  various  acceleration  increments  for 
cruise  altitude.  Fron  this  graph  the  exceedances  of  various  gust  velocities  can  be  read  off  using 
the  value  of  Uae/An  =  36ft/sec/g  given  by  Bullen  (1966),  and  converted  to  exceedances  per 
1000  nautical  miles  for  a  true  airspeed  of  270  kn.  The  results  are  given  in  Table  A2  at  the  end  of 
this  Appendix. 


A. 3  Bristol  Freighter  (SAFE,  New  Zealand) 

Reference:  Heath-Smith  (19S8) 

The  table  of  gust  exceedances  given  in  the  reference  was  normalised  (by  dividing  the  numbers 
of  exceedances  by  the  distance  flown)  and  a  simple  extrapolation  (assuming  the  logarithm  of 
the  number  of  exceedances  to  vary  linearly  with  the  gust  magnitude)  was  carried  out  to  obtain 
the  exceedances  of  the  innermost  levels.  The  result  is  given  as  Table  A3. 


A. 4  Super  Constellation  (Australia) 

A. 5  Super  Constellation  (Indian  Ocean) 

A. 6  Super  Constellation  (Far  East) 

A. 7  Super  Constellation  (Pacific  Ocean) 

Reference:  Heath-Smith  (1959) 

The  gust  exceedances  in  the  above  four  regions  have  been  tabulated  in  the  reference.  After 
division  by  distance  flown  they  have  been  listed  in  Tables  A4  to  A7  at  the  end  of  this  Appendix. 


A. 8  Viscount  (TAA) 

Reference:  Visick  (1965) 

The  reference  gives  a  table  of  gust  occurrences  in  4  ft/sec  gust  velocity  bands.  This  table 
was  integrated  and  interpolated  to  give  the  table  of  gust  exceedances  shown  as  Table  A8  at  the 
end  of  this  Appendix. 


A.9  Canberra  (TOPCAT) 

Reference:  Bullen  (1966) 

Bullen  (1966)  quoting  Wells  (1966),  gives  a  table  of  exceedances  of  gusts  for  routine  search 
flights  at  36,000  ft  in  the  TOPCAT  project.  The  relevant  results  are  repated  in  Table  A9. 


A.  10  Boeing  727  (Ansett) 

A. 11  Boeing  727  (TAA) 

Reference:  Hunter  and  Fetner(1967) 

Hunter  and  Fetner  show  a  table  of  gust  occurrences  in  4  ft/sec  bands  of  gust  velocity.  The 
calculation  of  gust  velocity  follows  the  American  practice  using  the  Pratt  and  Walker  gust 
alleviation  factor  and  an  old  formula  for  the  lift  curve  slope.  To  make  the  data  compatible  with 
the  other  data  presented  in  this  report  the  gust  velocities  were  factored  by  various  amounts, 
as  shown  in  the  table  below.  At  the  same  time  the  wing  lift  curve  slope  was  recalculated  using 
the  USAF  stability  and  control  DATCOM  procedure.  (NTIS  Accession  Number  N76-73204.) 
The  data  used  and  the  resultant  values  of  t/de/g  are  shown  below: 


W  =  57,652  kg 
5  =  148-5  mz 
c  =  4  •  32  m 


/R  =  7-5 

Leading  edge  sweepback  =  36° 
Mid-chord  sweepback  =  27° 

Section  lift  curve  slope  =  0-95  *2w 
Wing  thickness  =11%  chord 


mi 


Height 

(ft) 

CAS 

Hunter  and  Fetner 

British  method 

kn 

m/s 

a 

(per  rad) 

IWg 

(m/s) 

a 

(per  rad' 

Ut.lt 

(m/s) 

2,500 

185 

mm 

4-38 

19-48 

4-37 

18-49 

7,500 

304 

mSm 

4-67 

11-01 

4-67 

10-69 

12,500 

324 

167 

4-84 

9-8 

4-85 

9-62 

17,500 

332 

171 

5  02 

9-1 

5-04 

9-01 

22,500 

336 

173 

528 

8-46 

5-29 

8-48 

27,500 

324 

167 

5-48 

8-36 

5-48 

8-44 

32,500 

303 

156 

5-63 

8-60 

5-61 

8-75 

37,500 

282 

145 

5-82 

8-88 

4-37 

11-78 

In  Tact  the  two  calculations  of  gust  velocity  per  unit  acceleration  increment  give  very  similar 
values.  The  biggest  difference  (at  the  highest  altitude)  happened  mainly  because  the  Mach 
number  was  near  the  force  break  Mach  number  and  the  older  formula  for  the  lift  curve  slope 
did  not  allow  for  this.  After  the  necessary  corrections  and  interpolations  the  tables  of  gust 
velocity  exceedances  were  obtained  as  shown  in  Tables  AIO  and  Alt. 
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FIG.  1  QUANTITY  OF  DATA  AVAILABLE  AT  VARIOUS  ALTITUDES  FROM  THE 
VARIOUS  MEASUREMENT  PROGRAMMES. 

(The  area  of  each  rectangle  is  proportional  to  •  weighting  factor  (shown  in  Table  1) 
multiplied  by  the  distance  flown  in  the  height  band.) 
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FIG.  3  VARIATION  OF  lj0  WITH  ALTITUDE 

Australian  data  compared  with  ESDU  design  curve 


lIPGUSTS/OOWNCUSTS 


77 


# 


«> 


•* 


FIG.  5  VARIATION  OF  GUST  FREQUENCY  WITH  GUST  MAGNITUDE. 

Australian  data  compared  with  ESDU  design  curves  for  aircraft  without  cloud  warning  radar. 
All  Australian  measurements  in  height  range  are  pooled. 
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Fig.  6  VARIATION  OF  GUST  F  REQUENCY  WITH  GUST  MAGNITUDE. 

Australian  data  compared  with  ESDU  design  curves  for  aircraft  without  cloud  warning  radar. 
All  Australian  measurements  in  height  range  are  pooled. 
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